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We study RI/MOM renormalization constants of bilinear quark operators for Nf ~ 4 and the 
strong coupling constant for Nf = 2 + 1 + 1 using Wilson twisted-mass fermions. We use the 
"egalitarian" method to remove // (4) hypercubic artifacts non-perturbatively, which enables us to 
study physical quantities in a wide range of momenta. We then apply OPE in studying the running 
behavior of Zg and Us, from which we are able to extract the Landau gauge dimension-two gluon 
condensate (A^) which is of phenomenological interest. 
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1. Introduction 

The most systematic approach to study Quantum Chromodynamics (QCD) non-perturbatively 
is lattice QCD. By using the lattice formalism, one is often obliged to break some symmetries, 
among which is the continuum rotation symmetry. In Euclidean space, the 0(4) rotation symmetry 
is broken to H{A) or //(3) hypercubic symmetry depending on whether the lattice setup is the same 
on spatial and temporal directions. As a consequence, there are lattice artifacts which are only 
H{A) invariant but not 0(4) invariant. This is particularly an issue for computations of quantities 
hke renormalization constants since the associated statistical errors are often quite small, and the 
uncertainties from lattice artifacts become visible, thus deserve careful treatments. A popular solu- 
tion is to use the "democratic cuts" to select data points with relatively small H{A) lattice artifacts. 
Another approach, which we call the "egalitarian method" is to include the lattice artifacts 
explicitly in the data analysis. This approach allows one to use a wider range of data points and 
extract information from the lattice simulations more efficiently. In this work, we apply the "egali- 
tarian method" to calculations of renormalization constants of bilinear quark operators in RI/MOM 
scheme and the strong coupling constant in MOM-type Taylor scheme, obtaining results in a wide 
range of momenta with lattice artifacts under control. 

On the other hand, in the low and intermediate momenta region, the perturbative running for- 
mula may not be applicable and in principle one should take into account non-perturbative effects. 
It is suggested by Operator Product Expansion (OPE) that there could be various contributions from 
condensates. In our case, as we are working in Landau gauge, we claim that there could be con- 
tributions from the dimension two gluon condensate (A^) to some renormalization constants and 
running coupling constant a^. This quantity itself, (A^) in Landau gauge, is of great phenomeno- 
logical interests and has been studied extensively in literature (see and references therein). Here 
we show two independent methods to determine its value on the lattice. In both cases our results 
show strong evidences for a positive value of (A^) in Landau gauge. 

2. RI/MOM renormalization constants for Nf = 4 

In general, renormalization constants are needed to convert bare quantities computed on the 
lattice to the renormalized ones. They are essential for giving meaningful and accurate physical 
results. Renormalization constants can be calculated either perturbatively or non-perturbatively on 
the lattice. The perturbative approach suffers from systematic errors from truncations in perturba- 
tion series which are hard to estimate, while non-perturbative approaches suffer from discretization 
errors which could be systematically improved in practice, so generally the non-perturbative ones 
are preferred. In this work we adopt one particular non-perturbation method, the "RI'/MOM" 
scheme We focus on bilinear quark operators in the form & = where the Dirac matrix F 
represents /, 75 , , 7^1 75 , 7Ai 7v(^> v) ■ 

The RI/MOM scheme is defined by requiring the amputated green function Aff{p) equal to its 
tree level value at certain kinematic point: 

Z^'Z^^Tr[AM^TiPrX^-=,^- = h Z,(At' = p') = T^Tr[5,;,, (;.);/]. (2.1) 
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Here Zq is the quark wave function renormalization constant and Zff is the renormalization constant 
for operator ff. The amputated green function Kff{p) is defined by 

S{p)=a''Y^e-'P^{^lf{xmO)), Ke{p) = S-\p)Ge{p)S-\p) (2.2) 

X 

G{p) = a^l^e-'P^''-yH^if{xmO)^^ifm{y)). (2-3) 

x,y 

where the basic ingredients are quark propagators S{p) and Green functions G{p). 
On the lattice, four-momenta are discretized, i.e., they take values 

liiiK Nl Nl . ^ ^ ^ (2?i4 + l)7r A^r A^r 

Pi = in— «i = — ' = 1,2,3 p4 = n4 = — . (2.4) 

Nia 2 2 Nto 2 2 

Note that we have used anti-periodic boundary conditions for fermions in the temporal direction. 
First we compute quark propagators and vertex functions for a wide range of momenta = 
'X'"4~~T"''" 'T"- Then we apply the "egalitarian method" treatment to all data points 
to remove the hypercubic artifacts. The basic idea of this approach is that any polynomial function 
of p which is invariant under H{3) transformation is a polynomial function of the H{3) invari- 
ants /7^3 = L^^iP^i " = 2,4,6. Here we are using H{3) because we use different lengths and 
boundary conditions on temporal and spatial directions. Alternatively, one can use the functions of 
H{3) invariants as functions of H4 invariants p^l^\ = /?^, n = 2,4,6,8 plus the H{4) to H{3) 
symmetry breaking term oc {^p^ — p^j^). In practice, we find that the H{4) to H{3) breaking term is 
always small and can be safely ignored in the analysis. We thus use the ansatz that renormalization 
constants can be expressed as functions of H{4) invariants 

„[4] 

z'""{ap,a^KlcD) = Z^y^-'^'^a^ p^ap^^a^hlcD) +R{a^ p\a'^lcD)a^^ + • • • , (2.5) 

where Z'"" is the raw data on the lattice and z'^yp-'^^^r jj^g hypercubic corrected quantity with 
(almost) no H{4) lattice artifacts. Note that there could still be 0(4) invariant lattice artifacts in the 
quantity z'^yp-""''' and we will deal with them later. 

Gauge configurations are fixed to Landau gauge for calculating quark propagators. Under 
this gauge, OPE suggests that the dominant non-perturbative contribution to comes from the 
vacuum expectation value of (A^) = {A^^A"^). This term should be taken into account when we run 
Zq defined at different renormalization scales to some high scale, e.g., lOGeV, where perturbation 
theory is applicable for converting from RI/MOM scheme to MS scheme. Basically, we combine 
the perturbative running with OPE expansion in a single formula 

z^yp---ip^)=zp-'{p\^^) + +^.2p2«V, (2.6) 

where Z^"'' is the perturbative running which is available up to four-loop order C^ji^„^ is the 
Wilson coefficient for (A^) which is available up to ^(oi^) in certain schemes p]. The last term in 



Eq. (2.6) is the remnant discretization artifact which is invariant under 0(4) symmetry and is left 



over from the hypercubic corrections. 
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2.1 Lattice setup 

Our goal is to calculate the renormalization constants for Nf = 2 + 1 + 1 twisted-mass 
ensembles generated by the European Twisted Mass (ETM) Collaboration. As we are working in 
the mass-independent scheme, in which all the renormalization constants need to be extrapolated 
to the chiral limit, the strange and charm quarks used in these configurations are too heavy for our 
purposes. Instead, ETMC has generated dedicated Nf = 4 ensembles with four degenerate (twisted- 
mass) light quarks. It turns out that for Nf = 4 ensembles significant amount of work is needed to 
tune the action to the maximal twisted point |0] for achieving automatic i^(a) improvement. In 
practice, another solution was used: simulating with a pair of ensembles with equal absolute values 
of mpcAC (or 6) but opposite signs, then removing the &{a) artifacts by averaging the quantities 
obtained from these two ensembles. In this way, one avoids the work of fine tuning with the price 
of doubling the number of ensembles, and the strategy is indeed feasible as shown by David in 
their preliminary studies. 

To summarize, the procedures we need to go through for a complete analysis of the renormal- 
ization constants consist the following: 

• 0-average to remove &{a) artifacts. 

• Valence chiral limit to remove possible Goldstone-pole contributions for some RCs. 

• "Egalitarian" hypercubic corrections to remove discretization artifacts, i.e., &{a^^) 
terms. 

• Perturbative running plus OPE to run the renormalization constants to some high energy 
scale, e.g., lOGeV. The 0(4) invariant lattice artifacts are removed in this step. 

• Sea quark chiral Umit to remove remnant -dependence [^. 

Since the spacing is limited, here we only show one example of calculations of Zq, leaving a 
complete analysis of all RCs to later publications. We use one pair of 24^ x 48 ensembles with j8 = 



1.95 corresponding to lattice spacing a ^ 0.08 fm. In Fig. ( |l(a)D we show the so-called "fishbone" 
structure for raw data of Zq and hypercubic corrections with the "egalitarian" treatment. It can be 
clearly seen from the plot that even with the selection of "democratic" points the H{4) artifacts 
in the raw data are still significant (~ 5%) and could give misleading results without hypercubic 
treatments. In Fig. ( 1(b) ) we show the running fit of Zq{p^) to formula in Eq. (2^) where we have 
included the contributions from (A^) and 0(4) invariant lattice artifacts. 

After converting to physical units, we obtain the values of Z^""" and the gluon condensate (A?) 
at the scale of lOGeV: 

(At = lOGeV) = 0.7423 ±0.0026, (2.7) 
g^(A^)cMU=ioGev =4.44±0.12GeV2. (2.8) 

where only statistical errors are quoted. 

3. Running coupling constant for Nf = 2 + 1 + 1 

The parameter AgcD is a fundamental quantity of QCD and it is studied extensively in litera- 
ture. Its value can be determined by studying the renormalized running coupling constant from 



4 



RI/MOM RCs and the strong coupling constant 



X. Du 



.9=<ratb i — i — i 
.8=: ratio <. 9 - 
.7=<ralio<S ; - • -: 
.6 =< ralkj < .7 
.5 =< ratb < .6 
.4 =c ratio < 5 
.3=cratio<4 -■■ ■■ 

.2=<ratio<^ 

hvD corrected OWF 




(a) Hypercubic corrections 



(b) Running 



Figure 1: The left plot shows hypercubic corrections of Z^. Different colors represent data points with 
different values of ratio = p'^' / (p^)^- The yellow data points are those with relatively small H{4) artifacts 



and often selected by the "democratic" method. The right plot shows the running of Zq according to Eq. ( 2.6 ) 



experiment or from lattice calculations. Among many methods to calculate on the lattice, one 
category of approaches is to study the momentum behavior of Green functions. Particularly, it was 
shown in Ref. [§| that one way to study the running coupling constant is to study the ghost-ghost- 
gluon vertex in Landau gauge in the so-called Taylor scheme: 



471 A^oo 471 



(3.1) 



where A^) and G(/x^, A^) are the ghost and gluon dressing functions respectively. The advan- 
tage of this approach is that only two-point functions are involved thus one expects a good signal 
to noise ratio. This method has been successfully applied in previous Nf = iQ] and Nf = 2 [||] 
studies of Aqcd- Here we apply the same method to Nf = 2 + 1 + 1 ensembles where we include 
the dynamical strange and charm quarks as well as two degenerate light quarks. 

3.1 Lattice setup 

We use four Nf = 2 + 1 + 1 ensembles generated by ETMC, with different choices of j8, quark 
masses and volumes. The relevant parameters are listed in Table. (|T]). For generating these ensem- 
bles, Iwasaki action is used for the gauge part and Wilson twisted-mass action is used for both the 
degenerate light quark doublets and heavy quark doublets. The quark action is tuned to maximal 
twist by finding the value of Kent for which the PCAC mass vanishes mpcAC = 0- Consequently the 
physical (parity-even) quantities are automatic ^(a) improved [10]. 





'Q;rit 


a}Xi 


alia 


alls 


{L/af xT/a 


confs. 


1.95 


0.1612400 


0.0035 


0.135 


0.170 


32^ X 64 


50 




0.1612400 


0.0035 






48^ X 96 


40 




0.1612360 


0.0055 






32^ X 64 


50 


2.1 


0.1563570 


0.0020 


0.120 


0.1385 


48^ X 96 


40 



Table 1: Four ensembles used in our analysis. The one with j3 = 1.95 and volume 48^ x 96 is only used to 
check finite- volume effects. 
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(a) Hypercubic corrections 



. [i=l.95. a^|=0.0055 

. P=2.l.a(i,=0 0020 

■ p=l .95. a(i|=0.0035 

— OPE 

-■ g2<A>=0 




a(1.95)p 



(b) Running 



Figure 2: The left plot shows hypercubic corrections of the ghost dressing function. The right plot shows 



the running of aj according to Eq. (3.2) 



On each ensemble, we compute ghost and gluon dressing functions for a large range of mo- 
menta and apply the "egalitarian" hypercubic corrections to remove H{A) artifacts. This step, 



shown in Fig. ( 2(a) ) for the ghost dressing function, is essential since the perturbative logarithm 
scale dependence would be missing in the large momentum region {a^p^ > 1) had we used the data 
from "democratic" selections. 

After hypercubic corrections, we then calculate the running coupling constant in Taylor 



scheme using Eq. (3.1). It can be shown that with careful calibrations of lattice spacings [11], the 
data obtained from different ensembles superimpose with each other if we express in terms of 
momenta in units of one particular lattice spacing (a^^i 95 is used here). 

Again with the argument from OPE, we claim that there could be non-perturbative contribu- 



tions from a landau gauge {A ) to aj. Analogous to Eq. (p.q), we fit the running coupling constant 



at different scales to the formula [11] 



a 



pert 





1) 



(3.2) 



where aj^'^ is available to four-loop level in perturbation theory 
be calculated using the Wilson coefficient which is available to i?{a'^) 



The coefficient R{n^,ql) can 



The running fit for is shown in Fig. ( [2(b)[ ). The solid line represents the fitted curve to the 
running formula Eq. (O). The dotted line represents the value of aril^^) from pure perturbative 



running if we use the PDG value a{Z^) as input and run to lower scales. One can clearly see the 
discrepancy between the pure perturbative running and the running including OPE contributions. 
Furthermore, a careful analysis [jTl|] show that the difference between the lattice data and the fitted 
curve to Eq. (3.2) is scaling with 1 / p^, which seems to indicate that the next-to-leading correction 
in OPE is dominated by a /^{l/p^) term. This is certainly a non-trivial result and needs further 
investigations. 

By fitting to Eq. (^^), we obtain Aqcd and g^{A^) in the Taylor scheme. Converting to MS, 
we have the following results for j3 = 1.95 lattices: 



A^^ = 316± 13MeV, g^A^)R,^=ioG.Y = 4.5 ±OAGeY\ 



(3.3) 
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For the strong coupling constant, we can use the perturbative formula to run it to the scale of 
M20 with matching at the bottom quark mass threshold. The result is Oj^^iyijo) = 0.1 198(9) (10), 
where the first error is statistical and the second error is the estimated uncertainty from light quark 
mass effects. This result agrees well with the PDG value ^[^(Mzo) = 0.1 184(7) [ [l2| within errors. 



More details of our analysis can be found in Ref . [11] 



4. Conclusions and outlook 

In this work we studied the RI/MOM renormalization constants for bilinear quark operators at 
Nf = 4 and the running coupling constant at Nf = 2 + 1 + 1 using Wilson twisted-mass fermions. 
We have demonstrated that the "egalitarian" method is an effective way to eliminate hypercubic 
artifacts and its application is essential for obtaining reliable data in a wide range of momenta. 
Furthermore, we show that OPE plays an important role in the momentum range considered in our 
analysis. By taking into account the non-perturbative 0'{\/p^) contributions, we are able to extract 
the values of the Landau gauge gluon condensate iA^) from two independent approaches. In both 
Nf = 4 and Nf = 2 + 1 + 1 cases, we obtain values of g^(A^)s.^=ioGeV close to 4 or 5GeV^, which 
strongly suggest the existence of a non-zero Landau gauge dimension-two gluon condensate. The 
next step is to finish the analysis of renormalization constants for other bilinear and possibly twist-2 
quark operators. It would also be interesting to study the strong coupling constant in Nf = 4 so that 
one can investigate the dependence of (A^) on the number of sea quarks. 
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